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^ \ ABSTRACT 

0^' Wc demonstrate that the large scatter in the ultra-violet (UV) colours of intermediate- 

mass early-type galaxies in the local Universe and the inferred low-level recent star 
formation in these objects can be reproduced by minor mergers in the standard ACDM 
cosmology. Numerical simulations of mergers with mass ratios < 1:4, with reasonable 
^ \ assumptions for the ages, metallicities and dust properties of the merger progenitors, 

produce good agreement to the observed UV colours of the early-type population, if the 
infalling satellites are assumed to have (cold) gas fractions > 20%. Early- types that sat- 
isfy {NUV ~ r) < 3.8 are likely to have experienced mergers with mass ratios between 
1:4 and 1:6 within the last ^ 1.5 Gyrs, while those that satisfy 3.8 < {NUV — r) < 5.5 
' are consistent with either recent mergers with mass ratios < 1:6 or mergers with higher 

, mass ratios that occurred more than ~ 1.5 Gyrs in the past. We demonstrate that the 

early-type colour-magnitude relations and colour distributions in both the UV and 
optical spectral ranges are consistent with the expected frequency of minor merging 
activity in the standard ACDM cosmology at low redshift. We present a strong plausi- 
' bility argument for minor mergers to be the principal mechanism behind the large UV 

I scatter and associated low-level recent star formation observed in early-type galaxies 

in the nearby Universe. 
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1 INTRODUCTION 

The formation histories of early-type galaxies remains a con- 
troversial and often-debated topic in modern astrophysics. 
The bulk of the past effort on early-type systems has fo- 
cussed, almost exclusively, on their optical properties. The 
optical colours of the early-type population are predomi- 
nantly red, implying that the bulk o f the stellar mass in these 
systems forms at high redshift (e.g. iBower et al.lll993 ) . Fur- 
thermore, high a-enhancement ratios in their stellar spectra 
indicate that this star formatio n plausibly takes pla ce over 
short (< 1 Gyr) timescales (e.g. iThomas et al.lll999l ). How- 
ever, a drawback of optical data is its lack of sensitivity to 
moderate amounts of recent star formation (RSF), within 
the last Gyr or so. The optical spectrum remains largely un- 
affected by the minority of stellar mass that forms in these 
systems at low and intermediate redshift {z < 1), which 
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makes it difficult to measure early-type star formation his- 
tories (SFHs) over the last 8 billion years with significant 
accuracy using optical colours alone. 

A first step towards probing early-type SFHs over this 
period is to quantify their RSF at z ~ 0. Rest-frame UV 
photometry provides an attractive photometric indicator 
of RSF. While its impact on the optical spectrum is rela- 
tively weak (and virtually undetectable, given typical obser- 
vational and theoretical uncertainties), a small mass frac- 
tion (< 3%) of young (< 1 Gyr old) stars strongl y affects 
the r est-frame UV spectrum shortward of 3000A (|Kavirail 
l2008h . Furthermore, the UV rer nains largely un affected by 
the age-metallicity degeneracy (jWorthevI Il994l ) that typi- 
cally plagues optical analyses (Kaviraj et al. 2007a), making 
it an ideal photometric indicator of RSF. 

In Figure [1] we demonstrate the sensitivity of the UV to 
small mass fractions of young stellar populations. We con- 
struct a model in which an old (10 Gyr old) population 
contributes 99% of the stellar mass (shown in red), with a 
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Figure 1. The sensitivity of tlie UV to small mass fractions of 
young stellar populations. We construct a model in which an 
old (10 Gyr old) population contributes 99% of the stellar mass 
(shown in red), with a 1% contribution from stars that arc 0.3 
Gyrs old (show in blue). The combined spectral energy distribu- 
tion (SED) is shown in black. The UV output of the combined 
SED comes purely from the young population (blue) and that 
the shape of the optical spectrum - which determines the optical 
colours - changes only very slightly. We also indicate the (g — r) 
and (NUV — r) colours of the combined population (black) with 
those of the purely old population (red). While the {g — r) colour 
changes by ~ 0.1 mag from that of a purely old population, the 
{NUV — r) colour changes by ~ 2.5 mags! 



Figure 2. TOP LEFT: The (NUV - r) colour-magnitude rela- 
tion of nearby (0.05 < z < 0.06) early-type galaxies drawn from 
the SDSS DR4. Note that the GALEX NUV filter is centred at ~ 
2300 A. TOP RIGHT: The {NUV -r) colours of S PSS early-types 
plotte d against their stellar masses, taken from iGallazzi et al.l 
l|2005l ). Note that the Gallazzi et al. masses have been released 
through the public Garching SDSS catalog which can be found 
here: http://www.mpa-garching.mpg.de/SDSS/DR4/ The nom- 
inal stellar mass adopted for the spheroid in our simulations 
(4 X lO"*^" Mq, see Section 2 below) is indicated using the solid 
line. BOTTOM LEFT: The (g - r) colour-magnitude relation of 
nearby (0.05 < z < 0.06) early-type galaxies drawn from the 
SDSS DR4. 



1% contribution from stars that are 0.3 Gyrs old (show in 
blue). The combined spectral energy distribution (SED) is 
shown in black. The UV output of the combined SED comes 
purely from the young population (blue) and the shape of 
the optical spectrum - which determines the optical colours 
- changes only very slightly. We also compare the {g — r) and 
(NUV — r) colours of the combined population (black) and 
those of the purely old population (red). While the {g — r) 
colour changes by ~ 0.1 mag from that of a purely old pop- 
ulation, the {NUV — r) colour changes by ~ 2.5 mags! 

A new generation of early-type studies that have ex- 
ploited UV data from the GALEX space telescope (Martin 
et al. 2005) have shown that, in contrast to their optical 
colour-magnitude relations (OMRs), nearby (0 < z < 0.11), 
luminous {M{V) < —21) early-types show a large spread 
in their UV colour distribution of almost 5 mags - a direct 
consequence of the sensitivity of the UV to small amounts 
of recent star formation that is d emonstra t ed in Figure [1] 
above. Following the early work of lYi et al] (|2005l ). Kaviraj 
et al. (2007b; K07 hereafter) have demonstrated that, while 
a negligible fraction of the early-type population has pho- 
tometry consistent with no star formation within the last 2 
Gyrs, at least 30% show unambiguous signs of RSF, with 
stellar mass fractions of 1-3% forming within the last Gyr, 
with luminosity-weighted ages of ~ 300 — 500 Myrs. 

In the top left-hand panel of Figure [2] we show the 
{NUV - r) CMR of nearby (0.05 < z < 0.06) early-type 
galaxies drawn from the SDSS DR4. Note that the GALEX 
NUV filter is centred at ~ 2300 A. The large spread in 
the UV colours is in contrast to the small spread (a few 
tenths of a mag) in the optical {g — r) colour (bottom left- 
hand panel). Following K07, early-type morphology is es- 



tablished by visually inspecting each individual object us- 
ing its SDSS image. The galaxies shown in this figure are 
selected to have r < 16.8, since this is the redshift range 
within which visual inspection can be robustly performed 
using SDSS images (see Section 2 in K07). Note that the lu- 
minosity of a typical early-type galaxy, in the SDSS r-band, 
lies around M(r)* ~ —21.15 (see Figure 2 in Bernardi et al. 
2003). In the right-hand panel of this figure we present the 
{NUV — r) colours of the SDSS earl y-types pl otted against 
their stellar masses, taken from iGallazzi et al. I (|2005 ). Note 
that the Gallazzi et al. masses have been released through 
the public Garching SDSS catalog which can be found here: 
|http://www.mpa-garching. mpg.de/SDSS/DR4/ The nomi- 
nal stellar mass adopted for the spheroid in our simulations 
is 4 X lO'^" M0 (see also Section 2 below) where the scatter in 
the UV CMR becomes broad. This value is indicated using 
a dotted line. 

Note that, while the effects of RSF will be present in 
all diagnostics of star formation (including e.g. the optical 
{g — r) colour), the sensitivity of a particular indicator de- 
pends on the proportional change due to RSF compared to 
the typical uncertainties in that indicator. These uncertain- 
ties come both from measurement errors and from uncer- 
tainties in stellar models that are used to convert spectro- 
photometric quantities (e.g. colours) into estimates of phys- 
ical parameters such as the ages and mass fractions of young 
stars. Typical model and observational errors in optical fil- 
ters are ~ 0.05 mags (Sukyoung Yi priv. comm., see also 
Yi 2003) and at least 0.02 mags (including calibration un- 
certainties) respectively. The total resultant uncertainty in 
the optical {g — r) colours is ~ 0.08 mags (compared to a 
spread in this colour of ~ 0.3 mags from Figure [2]). In a 
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similar vein, the typical model and theoretical errors in the 
NUV magnitudes are both around 0.2 mags, yielding uncer- 
tainties in the {NUV — r) colours of ~ 0.3 mags (compared 
to a spread in this colour of almost 5 mags). The enhanced 
sensitivity of the NUV to RSF stems from the fact that the 
spread in the NUV colours is significantly (possibly an or- 
der of magnitude or more) larger than the uncertainties in 
the colours. While RSF does leave an imprint on the optical 
spectrum, the overall sensitivity of the [NUV — r) colour 
is much greater, making it a much more stringent test of 
the minor merger hypothesis than can be performed using 
optical colours alone. 

While the UV has revealed the unexpected presence of 
widespread low-level RSF in the local early-type population, 
past efforts have only measured the star formation activity 
without exploring the physical mechanism for RSF in these 
galaxies. In this paper, we explore the potential role of minor 
mergers in reproducing the UV colour-magnitude relation of 
intermediate-mass early-type galaxies at z ~ 0. It is worth 
noting that, although low level RSF can produce blue UV 
colours, only stars formed in the last ~Gyr contribute signif- 
icantly to the UV flux. Hence, an event like a minor merger 
- where one expects an age profile in the recent star for- 
mation - may not automatically produce blue UV colours, 
even if the supply of cold gas is adequately high to produce 
the observed mass fraction of young stars. In this study, we 
combine the results of numerical simulations of minor merg- 
ers with standard stellar models to probe the photometric 
properties of such events and test consistency with the ob- 
servations, given reasonable assumptions for the properties 
of the merger progenitors and the predicted frequency of 
merging activity in the ACDM paradigm. 



2 SIMULATIONS 

In this section we describe the numerical methodology used 
to study minor mergers between a typical elliptical galaxy 
and a satellite, where the mass ratio of the system is between 
1:4 and 1:10. A more complete description will be provided 
in Peirani et al. (in prep). Note that the nominal stellar mass 
of the spheroid used to describe the photometric properties 
of the simulations in Section 3 matches the region of the 
{NUV — r) vs. mass plot (see Figure [2] above) where the 
scatter is broadest. 



2.1 Initial conditions 

The elliptical is modelled using spherical dark matter (DM) 
halo and s tellar components, with a total mass of IO^^Mq. 
Following Ijiang fc KochaneM l|2007h . stars contribute 4% 
of this value. T he DM halo follows a Hernquist profile 
l|Hernauistlll990l ). with parameters cho sen so that the inne r 
part coincides with an NFW profile (|Navarro et al.l Il996l ) 
with a virial radius r2oo = 206 kpc and a concentration 
parameter C = 10, as indicated by p revious cosmological 
N-body simulations (|Dolag et al.ll2004 ). A Hernquist profile 
reproduces the 'de Vaucouleurs' 7?^''* surface brightness pro- 
files of typical elliptical galax;ies. The effective radius of the 
projected brightness is Ve = 4.3 kpc. 

The satellite is constructed using a spherical DM halo 
(with a Hernquist profile) which contains a disk containing 



stars and gas but no bulge. The mass of the disk repre- 
sents 5% of the total mass, with a gas fraction of either 
20% or 40%. The gas fraction s are consistent w ith the ob- 
served values from th e SDSS (lKannappan|l2004l ). Satellites 
are created following ISpringel et al.l 1 20051 ) and their rota- 
tion curves satisfy the baryonic TuUy-Fisher relation. In 
all si mulations satellites are pu t on prograde parabolic or- 
bits (jKhochfar fc BurkertllioO^ . with a pericentric distance 
Rp = 8 kpc, and initial separations of 100 kpc. 



2.2 Numerical method 

The s imulations are performed using the public GADGET2 
code jSpringelllioOsI ') with added prescriptions for cooling, 
star formation and feedback from Type la and II super- 
novae (SN). Approximately 380,000 particles are used for 
each experiment and the particle masses and gravitational 
softening lengths (e) involved are summarized in the table 
below. 





DM 


gas 


star (disk) 


star (E) 


Mass (lO^A^o) 


30.3 


4.5 


4.5 


13.5 


e (kpc) 


0.4 


0.5 


0.5 


0.4 



The cooling and star formation (SF) r ecipes fol- 
low the pr e script ions of iThomas fc CouchmanI (| 19921 ') and 
iKatz et al] (| 19961 ') respectively. Gas particles with T > lO^'K 
cool at constant density (with the assumption of solar metal- 
licity) for the duration of each timestep. Gas particles with 
T < 2 X 10* J\, number density n > 0.1 cm~'^, overden- 
sity Spgas > 100 and V.w < form stars according to the 
standard star formation prescription: dpt/dt = Ctpgas/tdyn, 
where p* refers to the stellar density, tdyn is the dynamical 
timescale of the gas and c, is the SF efficiency. Instead of 
creating new (lighter) star particles, we implement the SF 
prescription in a probabilistic fashion. Assuming a constant 
dynamical time across the timestep, the fractional change 
in stellar density, Ap*/p« = 1 — exp(— c* At/tdyn). For each 
gas particle, we draw a random number (r) between and 
1 and convert it to a star if r < Ap«/p«. 

The energy injection into the inter-stellar medium 
(ISM) from SN, which regulates the star formation rate 
(SFR), is modelled following the approach of Durier fc de 
Freitas Pacheco (2007, in prep.). Instead of assuming 'in- 
stantaneous' energy injection, we include the effective life- 
time of SN progenitors using the rate of energy injection 
HsN- For this, we consider stellar lifetimes in the mass 
ranges 0.8 Mq < m < 8.0 M© and 8.0 M© < m < 80.0 M© 
for Type la and Type II progenitors respectively. Using a 
Salpeter initial mass function for Type II SN gives: 



HsN,, = 2.5 X l0-^»(-^)£siv( 



1300 



r(Myr) 



erg.s 



(1) 

where Esn = lO^'^ erg, m, is the mass of the star particle 
and 3.53 < r < 29 Myr. For Type la SN, the heating is 
delayed, since they appear tp = 0.8— 1.0 Gyr aft e r the onset 
of star formation. Following Ide Freitas Pachecd (|l998t ). the 
probability of one event in a timescale r after the onset of 
star formation is given by: 
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Figure 3. The star formation rate of a minor major with mass 
ratio 1;6, where the gas fraction of the satellite is 20% and c» = 
0.05. The dashed line represents the same simulation with 10 
times more particles. At the 'final plunge' the satellite disappears 
into the parent elliptical, so that images (e.g. through an SDSS 
r-band filter) would indicate a single spheroidal object. 



HsM,^ = 4.8 X 10-^°(g)iJ^^(^)''%rg.s- (2) 

Eqns (1) and (2) are used to compute the energy re- 
leased {Ei) by SN derived from a star particle i, and a frac- 
tion 7 of this energy is deposited in the j"" neighbour gas 
particle by applying a radial kick to its velocity with a mag- 
nitude Awj = ^ (2wj'yEi/mj), where Wj is the weighting 
based on the smoothing kernel and rrij is the mass of gas 
particle j. We note that all gas neighbours are located in 
a sphere of radius Rsn, centered on the SN progenitor, to 
avoid spurious injection of energy outside the SN's region of 
influence. 

In the simulations presented in this study, we use 7 — 
0.1, Rsn = 0.4 kpc and vary c« in the range 0.01 < c« < 0.1. 
When satellites are isolated, these parameters lead to a 
quasi-constant SFR between 0.2 and 3.OM0yr~^, depending 
on the mass and initial gas fraction of the satellite. These 
SFRs are in good agreement with those for low-mass ob- 
jects found in prev ious simulations of isolated galaxies (e.g. 
IStinson et al]|2006h . 

Finally, we have checked that increasing the resolution 
of the simulations (by a factor of 10) does not affect the star 
formation history or the derived colours so that the conclu- 
sions presented in this work are robust. A typical example 
of the SFH is shown in Figure [3] 



3 THE CASE FOR MINOR MERGERS 

We begin by exploring the predicted photometry from an 
ensemble of minor merger simulations (described in Section 
2) with the assumptions that (a) both merger progenitors 
have solar metallicity, (b) the elliptical and satellite have 
luminosity-weighted ages of 9 and 5 Gyrs respectively and 
(c) the dust extinction due to the ISM in the remnant is 
given by E'J% ~ 0.05. 

We briefly describe the motivation for modelling the 
underlying population of the spheroid using a 9 Gyr old 



simple stellar population (SSP). It is well established that 
the bulk of the stellar populations in elliptical galaxies are 
uniformly old. This is demonstrated by the fact that (a) they 
show red optical colours w ith small scatter (e.g. lBower et al] 
ll992l : IStanford et al.|[l99i ) and (b) they exhibit high alpha- 
enhancement (e.g. [Mg/Fe]) ratios, which implies that the 
star formation took place on timescales shorter than the 
typical o nset timescales of Ty pe la supernovae i.e. < 1 Gyrs 
(see e.g. I Thomas et al.l Il999j l . Hence the underlying popu- 
lation can be approximated by an old SSP. The particular 
choice of 9 Gyrs is motivated by the fact that Bernardi et 
al. (2003), who recently studied the SDSS elliptical galaxy 
population in the nearby Universe, were able to fit their op- 
tical colour-magnitude relations with an SSP with an age 
of 9 Gyrs. It is worth noting tha t opti cal colour evolution 
virtually stops after ~ 6 Gyrs (|Yill2003l ). This means that a 
6 Gyr old SSP looks very similar to an older stellar popula- 
tion e.g. 9 Gyrs. In other words we could replace the 9 Gyr 
SSP with a 6 or 10 Gyr SSP and our results will not change. 
Note that similar techniques (i.e. using an old SSP to rep- 
resent the underlying stellar population of elli pticals) have 
been freque ntly used by previous studies, e.g. iTrager et al.l 
l|2000h and iFerreras fc Silkl l|2000l ). Finally, stellar popula- 
tions that are older than ~ 2 Gyrs do not contribute to the 
UV (see the bottom panel of Figure 7 in K07). As a result, 
the underlying population of an elliptical galaxy will not 
contribute to the UV at all. 

The estimate for the ISM extinction is an average value 
for local early-types, derived by K07 from parameter estima- 
tion using GALEX (UV) and SDSS (optical) photometry. 
Since the youngest stars are expected to reside in molecu- 
lar clouds ( MCs) which have short lifetimes of a few tens 
of Myrs (e.g iBlitz fc Shulll980l : iHartmann et aLlbOOll ) and 
dust extinction several times larger than that due to the 
ISM alone (e.g. ICharlot fc Falll [2000f ) . we explore MC life- 
times in the range 0-50 Myrs and MC extinctions in the 
range 0.05 < E^Sv < 0.5. Thus, when the synthetic pho- 
tometry from the simulations is constructed, stars with ages 
less than the MC lifetime in question are subject to the pre- 
scribed MC extincti on. Th e stellar models used in this study 
are described inlxl (|2003t ). 

In Figure |4] we present the synthetic photometry from 
various configurations for mergers with mass ratios of 1:10. 
The remnant is 'observed' at the point where the satellite 
finally disappears into the parent elliptical, so that images 
(e.g. through an SDSS r-band filter) would indicate a single 
spheroidal object. The colours shown are therefore the bluest 
possible for each scenario where the system appears to be 
one object. Star formation declines after this 'final plunge' 
(see Figure [Sjl and the remnant reddens, in the {NUV — r) 
colour, by ~ 0.8 mag/Gyr for mass ratios of 1:4 and 1:6 and 
by ~ 0.5 mag/Gyr for a mass ratio of 1:10. 

We find that, with reasonable MC properties - e.g. life- 
times > 30 Myrs (red and green colours) and high dust ex- 
tinctions (large symbols) - 1:10 mergers can reproduce most 
but not all of the scatter in the observed UV colours. In 
particular, the bluest UV colours {NUV — r < 4) cannot be 
accounted for by 1:10 mergers because the RSF induced is 
inadequate, even when the satellite has a high gas fraction 
(~ 40%). The observed early-type UV colour-magnitude re- 
lation (CMR; shown using filled grey circles in Figure |4]) is 
restricted to the redshift range 0.05 < z < 0.06 (to ensure 
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Figure 4. Predicted photometry from an ensemble of minor 
merger simulations with mass ratios of 1:10. Both progenitors 
are assumed to have solar metallicity with the elliptical and the 
satellite having ages of 9 and 5 Gyrs respectively. The ISM dust 
extinction in the remnant is given by Eg^J'y ~ 0.05. Symbol types 
represent merger configurations - the star formation efficiency (c*) 
and gas fraction (G) of the accreted satellites arc shown in the top 
left legend. MC extinctions arc shown using colours, while sym- 
bol sizes represent MC lifetimes. Grey dots represent the observed 
colours of luminous {L* or above) early-type galaxies in the red- 
shift range 0.05 < z < 0.06. Note that typical uncertainties in the 
observed colours are 0.2 mag and that the synthetic photometry 
has been redshifted to z = 0.065 for a direct comparison. 



Figure 5. TOP: Predicted photometry from an ensemble of mi- 
nor mergers with mass ratios of 1:10. We assume a spread in 
the age (1-9 Gyrs; shown using symbol sizes) and metallicity 
(O.2Z0 — 2Zq; shown using colours) of the satellite. The observed 
colours of the local, intermediate-mass (L* or above) early-type 
population is shown by the black dots. Large grey triangles in- 
dicate the positions of simple stellar populations with half-solar 
and solar metallicities that form at z = 3. 



robust morphological classification from SDSS images) and 
r < 16.2 (which corresponds to galaxies more luminous than 
L* at this redshift). 

Recalling that Figure |4] assumes solar metallicity and a 
single age for the satellite (5 Gyrs) , we now explore a wider 
parameter space where we vary the metallicity of the satel- 
lite in the range 0.2Zq to 2Zq and its luminosity-weighted 
age in the range 1-9 Gyrs. The analysis is restricted to real- 
istic MC properties - lifetimes > 30 Myrs and extinction > 

0. 3. We present these results in Figure[5l While the predicted 
photometry falls in the centre of the locus of observed galajc- 
ies in Figure 31 we find that a reasonable spread in the age 
and metallicity of the satellite can reproduce the 'horizontal' 
spread in the (u — r) colours. However, such a spread in age 
and metallicity cannot mimic low MC lifetimes/extinctions 

1. e. if we restrict our analysis to realistic MC properties only, 
low satellite ages and metallicities remain unable to repro- 
duce the bluest UV colours. We also indicate, in Figure [SI 
the early-type galaxy fractions bluer than {NUV — r) < 3.8 
(the colour limit of the 1:10 mergers with realistic MC prop- 
erties) and in the colour range 3.8 < NUV — r < 5.5. K07 
estimated that, within theoretical and observational uncer- 
tainties, galaxies bluer than {NUV — r) ~ 5.5 are very likely 
to have had some RSF (within the last ~Gyr), while galax- 
ies with {NUV — r) > 5.5 can be consistent with both RSF 
and purely old (> 2 Gyrs old) stellar populations. 

We now repeat the analysis with mergers that have mass 
ratios of 1:6 and 1:4. Figure [6] indicates that the blue end of 
the early- type UV colours, which is inconsistent with 1:10 



• 1 ;6 mergers 

• 1 :4 mergers 




• • » 



o o= 



(u-r) 



Figure 6. Predicted photometry from the same ensemble of mi- 
nor merger simulations as in Figure |4] but with mass ratios of 
1:6 (blue) and 1:4 (red). Both progenitors are assumed to have 
solar metallicity with the elliptical and the satellite having ages 
of 9 and 5 Gyrs respectively. The ISM dust extinction in the rem- 
nant is given by Eg^_^y ~ 0.05. We only show configurations with 
realistic MC properties - lifetimes > 30 Myrs and extinction > 
0.3 ^ 



mergers, can indeed be reproduced by 1:6 and 1:4 mergers 
with realistic assumptions for MC properties. 



6 Sugata Kaviraj et al. 



4 MERGER STATISTICS IN ACDM AND 
REPRODUCTION OF THE OBSERVED 
EARLY- TYPE COLOURS 

While the UV colours of intermediate-mass early-type galax- 
ies appear consistent with minor mergers, the reproducibil- 
ity of the entire UV OMR depends on the frequency of merg- 
ing activity at low redshift. In FigureQwe present the aver- 
age fraction of early-type galaxies (in dark matter halos of 
mass 10^^ M0 or above) that are predicted to have had one 
(solid), two (dashed), three (dot-dashed) and four (triple 
dotted) 1:X mergers in a given look-back time in ACDM. 
The value of 'X' is indicated in each panel. Merger trees 
are generated using the semi- analytical model of Kho chfar 
and Burkert (2003, 2005) and iKhochfar fc Silkl (|2006l ) and 
morphology is traced using stellar bulge:total (B/T) ratios 
- early-type galaxies are assumed to have B/T > 0.7. 

The merger fractions are calculated by dividing the 
number of mergers that occurred in look-back time (t) bins 
in the histories of spheroids by the number of such spheroids 
at 2; = 0. For example, the point at t = represents the 
merger fraction within look-back times of and 1 Gyr, while 
the point at t = 1 represents the merger fraction within 
look back times of 1 and 2 Gyrs and so on. The merger 
fraction increases with look-back time because the merger 
rate in the Universe incr e ases at higher redsh ift (see also 
e.g. iGottlober et al.ll200ll : iLe Fevre et~ai1l200Gl '). Note that 
the merger fraction increases monotonically until the merger 
rate in the Universe peaks and then drops off. 

Before comparing the photometric predictions from our 
model machinery to the observations, we briefly note that, 
while major mergers (which have mass ratios between 1:1 
and 1:3) might also be expected to contribute to the scatter 
in the UV CMR, a study of the GALEX NUV photom- 
etry of a samp l e of o ngoing major mergers, identified by 
iMcIntosh et al.l |2003) from the SDSS, indicates that ma- 
jor merger progenitors lie on the broad UV red sequence 
{NUV -r > 4.7). While a caveat here is that the Mcintosh 
et al. sample is small (compared to the SDSS galaxy popu- 
lation and especially after cross-matching with GALEX), it 
is reasonable to conclude that major mergers will contribute 
only to the broadness of UV 'red sequence' [NUV — r > 4.7), 
not to the full extent of the UV scatter and certainly not to 
the blue end of the UV CMR. 

We now check if the predicted LCDM minor merger ac- 
tivity (i.e. events with mass ratios between 1:4 and 1:10), 
convolved with the predictions from the basic set of numeri- 
cal simulations described above, is able to simultaneously 
reproduce the UV and optical CMRs of the low-redshift 
early-type population. Using the SDSS luminosity function 
of low-redshift early-type galaxies extracted in K07, we con- 
struct a large Monte-Carlo (MC) ensemble of 50,000 sim- 
ulated objects and compare the synthetic CMRs from this 
ensemble with the observed early-type CMRs in both the 
(NUV — r) and [g — r) colours. We assume (a) realistic 
MC properties i.e. random MC lifetimes between 30 and 50 
Myrs and 0.3 < E%'Sv < 0.5 (b) a gaussian distribution in 
the SSP-weighted satellite ages, which peaks at 5 Gyrs, with 
a width of 2 Gyrs and (c) a uniform distribution in the satel- 
lite metallicities in the range Q.2Zq — 2Zq. We only consider 
merger events within the last 2.5 Gyrs and note that includ- 
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Figure 7. The average fraction of spheroidal galaxies (in dark 
matter halos of mass 10^'^ Mq or above) that are predicted to 
have had one (solid), two (dashed), three (dot-dashed) and four 
(triple dotted) 1:X mergers in a given look-back time in ACDM. 
The value of 'X' is indicated in each panel and the dispersions in 
the fractions are ~ 20%. Fractions are calculated by dividing the 
number of mergers that occurred in look-back time {£) bins in the 
histories of spheroids by the number of such spheroids at 2; = 0. 
For example, the point at t = represents the merger fraction 
within a look-back time of and 1 Gyr, while the point at t = 1 
represents the merger fraction within a look back time of 1 and 
2 Gyrs and so on. The merger fraction increases with look-back 
time because the me rger rate in the Universe increases at highe r 
redshift (see also e.e. lGottlober et al.ll200ll:lLe Fevre et al.ll2000t) . 
Note that the merger fraction increases monotonically until the 
merger rate in the Universe peaks and then drops off. 



ing larger time windows does not alter our results because 
the UV flux decays rapidly after 2 Gyrs anyway. 

In Figure |8] we draw a random subset, equal to the num- 
ber of observed galaxies in our comparison redshift range 
(0.05 < 2 < 0.06), from the parent MC ensemble of 50,000 
objects and compare their [NUV — r) and (g — r) colours to 
that of the observed early-type population. A comparison 
between the colour histograms of the parent MC ensem- 
ble and the observed early-type population is also shown 
in the inset to each figure. Note that both histograms are 
normalised to 1. We find that, given the assumptions listed 
above, there is good quantitative agreement between the 
synthetic and observed CMRs and colour distributions in 
both UV and optical colours. In other words, the predicted 
minor merger activity in the standard model is able to repro- 
duce the UV/optical properties of the early-type population 
(in particular the distribution scatter to blue UV colours) 
with a satisfactory level of accuracy. This, in turn, provides 
a very strong plausibility argument for minor merging driv- 
ing the recent star formation observed in early-type galaxies 
at low redshift. 
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fewer pre-mergers being found in total. An added problem 
is that computing relative velocities requires spectroscopic 
redshifts. As a result smaller, fainter galaxies at the mag- 
nitude limit of typical spectroscopic surveys will not have 
redshifts and therefore minor mergers cannot be efflciently 
identified (since only the larger progenitor has a redshift). 
Hence the lack of large numbers of minor merger closepairs 
in typical observational studies is not a good indicator of 
the true frequency of such events in the nearby Universe. 

A better alternative is to look for post-mergers i.e. ob- 
jects that show tell-tale signatures of recent merger activ- 
ity. An added advantage of studying post-merger systems is 
that one is sure that the merger has actually taken place, 
whereas with pre-merger systems it is just a guess (albeit a 
well-educated one!). While a large observational survey like 
the SDSS could be the perfect test-bed for performing such 
a study, the typical 50 second exposure SDSS images are not 
deep enough to detect faint morphological disturbances from 
minor mergers. However, van Dokkum (2005) have recently 
used very deep optical photometry to show that over 70% of 
early-types on the optical red sequence show morphological 
signatures of merging. The morphological features, e.g. fans, 
tails, shells, are faint and red to surface brightness limits of 
^ ~ 28 mag arcsec"^) and consistent with minor mergers 
(where the induced star formation is at a very low- level). 
These features do not appear in the corresponding SDSS 
images and an analysis of the UV magnitudes of these red 
mergers (using GALEX) indicates that their UV CMR is 
similarly broad to what has been found for the general SDSS 
early-type population in K07 (Kaviraj and van Dokkum, in 
prep). 



Figure 8. TOP: Comparison of the simulated (NUV — r) colours, 
generated by convolving the photometric predictions from minor 
merger simulations with the predicted minor merger activity in 
LCDM (black), with the observed [NUV — r) colours of the low- 
redshift early-type population. Colour histograms (normalised to 
1) are shown in the inset. BOTTOM: The corresponding plot for 
the optical (g — r) colour. 



4.1 A note about observational signatures of 
merging 

It is worth noting here that the arguments presented above 
imply that an appreciable number of early-type systems in 
the nearby Universe should either be in a 'closepair' or 'pre- 
merger' system with a satellite or exhibit morphologies that 
are consistent with recent merging events. 

We should mention first that several observational stud- 
ies have looked at closep airs from a variety of surveys across 
a ran ge of redshifts (e.g. jPatton et al.|[200ol : iLe Fevre et al.l 
I2OOOD . While such studies are very useful, a potential com- 
plication is that the dete c tion o f closepairs relies on a range 
of criteria. iPatton et all (|2000l ). for example, tag pairs of 
galaxies as 'pre-mergers' where the tranverse separation on 
the sky is less than 20h~^ kpc and the relative velocities are 
less than 500 kms~^. The number of closepairs naturally de- 
pends on the criteria being used. For example, making the 
criteria stricter increases the likelihood that each candidate 
system is truly a 'pre-merger'. However it also results in 



5 DISCUSSION AND SUMMARY 

We have compared the UV colours of nearby (0.05 < z < 
0.06) early-type galaxies with synthetic photometry derived 
from numerical simulations of minor mergers, with reason- 
able assumptions for the ages, metallicities and dust prop- 
erties of the merger progenitors. Observation al estimates 
for sa tellite gas fractions have been taken from iKannappanI 
(l2004h and minor merger simulations have been performed 
using these gas fractions. We have then appealed to the 
merger statistics in the standard ACDM paradigm to check 
whether the minor merger activity could plausibly drive the 
scatter in the UV CMR at low redshift. 

We have found that the bluest end of the early-type UV 
CMR {NUV — r< 3.8) is consistent with mergers that have 
mass ratios between 1:4 and 1:6 (and cannot be reproduced 
by events with mass ratios less than or equal to 1:10), as- 
suming that the infalling satellites have gas fractions around 
^ 20% or higher, which are consisten t with the observa - 
tionally constrained gas fractions from IKannappanI (j2004[ ). 
Early-types with intermediate UV colours (3.8 < NUV —r < 
5.5) are consistent with either recent minor mergers with 
mass ratios less than 1:6 or mergers with higher mass ratios 
more than ~ 1 Gyr in the past. Major mergers are likely 
only to contribute to the broadness of the UV red sequence 
and not to the blue scatter in the UV CMR. Furthermore, 
we have demonstrated that the predicted minor merger ac- 
tivity in the standard model, convolved with photometric 
predictions from our fiducial set of numerical simulations, is 
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able to simultaneously reproduce the UV and optical CMRs 
and colour distributions of the low-redshift early-type pop- 
ulation, in particular the large scatter to blue colours. This, 
in turn, provides a strong plausibility argument for minor 
mergers being responsible for the RSF in early-type galax- 
ies in the nearby Universe. 

We note here that our study does not utilise a full-blown 
semi-analytical model for two important reasons. Firstly, the 
amount of gas available at late times in such a model de- 
pends on the baryonic recipes implemented within it. Dif- 
ferent models can predict different gas fractions in satellites 
at late times. Secondly, in current semi-analytical models 
the merger event is modelled 'instantaneously' i.e. the full 
age profile of stars formed in the merger is not taken into 
account. The correct way to implement minor merger events 
in a semi-analytical model is to apply star formation histo- 
ries from simulations such as those presented in this paper 
to merger events in the model. A full analysis with such 
an implementation will be provided in a forthcoming paper 
(Khochfar et al., in prep). In this paper we have appealed 
only to the merger fractions predicted by LCDM, which are 
robust. Finally, the primary reason for using the Khoch- 
far and Burkert model is that it uses the Extended Press 
Schechter (EPS) formalism to generate merger trees, result- 
ing in infinite mass resolution which is important because 
we are specifically looking at minor mergers involving small 
(satellite) objects. 

Finally, we note some possible caveats to the analysis 
presented here. Firstly, although it is reasonable to infer rela- 
tively high g a s frac tions (> 20%) for infalling satellites based 
iKannappanI (|2004l ). it is worth noting that the gas frac- 
tions are derived from calibrations calculated by combining 
SDSS photometry and Iff measurements. Confirmation of 
these gas fractions (especially in satellite galaxies) requires 
larger and deeper surveys that yield gas mass measurements, 
which are unavailable right now but may be possible us- 
ing future instruments such as Herschel. Secondly, we note 
that although cooling from hot gas halos could potentially 
contribute to the supply of cold gas that then drives star 
formation, this cha nnel is unlikely in galaxi es in the mass 
range studied here. iDekel fc Birnboiml l|2006l ) show that, at 
late epochs (z < 2), the gas in halos above a critical shock- 
heating mass (IO^'^Mq; consistent with the galaxies consid- 
ered in this study) is heated by a virial shock, leading to long 
cooling timescales that effectively shut down the gas sup- 
ply and subsequent star formation. Furthermore, unchecked 
accretion from hot gas halos wou ld result in early-ty pes be- 
coming too massive and too blue (|Benson et al.|[200a ) . Since 
there is no reason to believe that that z ~ is a preferential 
epoch for gas cooling, evidence from the luminosity function 
of massive galaxies renders star formation from cooling flows 
very unlikely. This, co mbined with t he abundance of minor 
merger features in the Ivan DokkumI (2005) study, strongly 
indicates that the recent star formation is driven not by gas 
cooling but by merging activity. 

While this study provides a strong plausibility argu- 
ment for minor mergers being the principal mechanism be- 
hind the large UV scatter and associated low-level recent 
star formation in early-type galaxies, similar studies are re- 
quired at intermediate redshifts (0.5 < z < 1) to check 
whether the evolution of this scatter is consistent with the 
ACDM merger statistics. In addition, deep/high-resolution 



imaging i.e. with the Hubble Space Telescope are required 
to confirm the coincidence of morphological signatures pro- 
duced by merging with UV excess in these galaxies. Current 
galaxy surveys such as COSMOS (Scoville et al. 2007) make 
such analyses possible and results from such studies will be 
presented in a forthcoming paper. 
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